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We calculate analytically properties of holographic superconductors in the probe limit. We ana- 
lyze the range 1/2 < A < 3, where A is the dimension of the operator that condenses. We obtain 
the critical temperature in terms of a solution to a certain eigenvalue problem. Near the critical 
temperature, we apply perturbation theory to determine the temperature dependence of the con- 
densate. In the low temperature limit we show that the condensate diverges as T~ A/3 for A < 3/2 
whereas it asymptotes to a constant value for which we provide analytic estimates for A > 3/2. We 
also obtain the frequency dependence of the conductivity by solving analytically the wave equation 
of electromagnetic perturbations. We show that the real part of the DC conductivity behaves as 
e -A g /T anc [ es timate the gap A 9 analytically. Our results are in good agreement with numerical 
results. 

PACS numbers: 11. 15. Ex, 11.25.Tq, 74.20.-z 



I. INTRODUCTION 

The BCS theory of superconductivity [l[ is a microscopic theory of weakly coupled superconductors and describes 
them with great accuracy. However in the last few decades it has been discovered that there is another class of 
superconductors that are strongly coupled. Conventional theories have failed to describe these systems. The pairing 
mechanism and the normal state of the system before the onset of superconductivity remain open questions (see [2| 
for an introduction to strongly coupled superconductors). In recent years it has been shown using the AdS^+i/CFTrf 
correspondence [|| that for d = 3 and d — 4 we can create a strongly coupled superconducting system from its 
gravitational dual 4-7], also including an external magnetic field [1,0. The d — 3 systems are especially interesting 
since many strongly coupled superconductors are effectively 2 + 1— dimensional with the conductivity occurring in 
planes. Phenomena such as the Hall effect [l(| and Nernst effect [lll - [l3j have also been shown to have dual gravitational 
descriptions. Other areas of Condensed Matter Physics have been discussed in terms of the AdS/CFT correspondence; 
sec hj for a review. The analysis for the most part has been done numerically with the exception of ref. [15| in which 
the WKB approximation was employed in the calculation of conductivity. 

Until recently the ground state of a holographic superconductor has been unreachable because a numerical solution 
to the non-linear field equations becomes increasingly cumbersome as the temperature approaches zero. Progress was 
made in [l6l [l7| where the ground state was reached numerically and the T — conductivity was analyzed. Using 
the ground state, the Fermi surface of these systems was then explored [18|. 

Here we explore the properties of holographic superconductors using analytical techniques aiming at elucidating, 
among other things, the nature of the ground state. We concentrate on the simplest case by working in the probe limit 
in which the back reaction to the bulk metric can be ignored. We obtain analytic approximations to the solutions of 
the non-linear field equations both near the critical temperature T c and in the low temperature limit T — ¥ 0. This 
enables us to avoid problems with numerical instabilities which one encounters when one solves the field equations 
numerically near zero temperature. 

The paper is organized as follows. In scction|TT]we review the field equations. In section Hill we discuss the properties 
of superconductors near the critical temperature. In section ITVl we analyze the zero temperature limit. In section fVl 
we discuss the conductivity. Finally, section IVT1 contains our concluding remarks. 
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II. FIELD EQUATIONS 



We are interested in the dynamics of a scalar field of mass m coupled to a U(l) vector potential in the backgound 
of a 3 + 1— dimensional AdS Schwarzschild black hole with planar horizon of metric 



ds 2 = -f(r)dt 2 



dr 2 



+ r 2 dx 2 



f(r)=r 2 - r -± 



(1) 



/(r) r 
in units in which the AdS radius is I = 1. The radius of the horizon is r + and the Hawking temperature is 

Assuming that the scalar field is of the form ty(r) and the potential is an electrostatic scalar potential the field 

equations are @ 



1 4 



Under the change of coordinates 
the field equations become 



T , / $ 2 TO 2 \ 
* + * 

V/ 2 fJ 

2 * 2 



■'■+ 



(3) 



(4) 



2 + z J 



$2 



>^(1-Z 3 ) 2 



2* 2 



z 2 (l-z 3 ) 



(5) 



where prime now denotes differentiation with respect to z, to be solved in the interval (0, 1), where z = \ is the 
horizon and z — is the boundary. 

Near the boundary (z — > 0), we have approximately 



(±).A± 



$ W /i ^-Z 



where 



A± = -± „ 
* 2 V 4 



(6) 



(7) 



While a linear combination of asymptotics is allowed by the field equations, it turns out that any such combination is 
unstable [2lj |. However, if the horizon has negative curvature, such linear combinations lead to stable configurations 
in certain cases [2^ . 

Thus, the system is labeled uniquely by the dimension A = A±. We shall examine the range 



i<A<3 



(8) 



where A > 3/2 (A < 3/2) if A = A + (A = A_), corresponding to masses in the range > m 2 > —9/4 (above the 
Breitenlohner- Friedman bound (l9l [2G|V 
Demanding at the horizon 



$(1) = 



(9) 



/j is interpreted as the chemical potential of the dual theory on the boundary, p is the charge density on the boundary 
and the leading coefficient in the expansion of the scalar yields vacuum expectation values of operators of dimension 
A±, 



(10) 



The field equations admit non-vanishing solutions for the scalar below a critical temperature T c where these operators 
condense. 
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III. NEAR THE CRITICAL TEMPERATURE 



At the critical temperature T c , & = 0, so the field equation ([3]) for the electrostatic potential reduces to <&" = 0. 
We may set 



*(*) = Ar +C (l - z) , \ = -%- 



(11) 



+c 



where r +c is the radius of the horizon at T — T c . 

As T — > T c , the field equation for the scalar field approaches the limit 



2 + z 3 



# = 



A 2 



z(l-z 3 )^ ' z 2 (l-z 3 )~ (1 + z + z 2 ) 2 



(12) 



which is valid even at T — T c , because it is linear in (one may divide "J by the leading coefficient in the expansion 
around z — before letting <f — > 0). 

To match the behavior at the boundary, define 



*(*) 



V2r$ ( ' 



where we used eq. (|10[) to express the leading order coefficient so that F is normalized to F(0) = 1. 
We deduce 



F 



2 + z 3 
1 - z 3 



A 2 



to be solved subject to the boundary condition 



1 , A 2 z 
- 2A F + F = - 

J 1-z 3 (1 + z + z 2 ) 2 

F'(0) = 



F 



The eigenvalue A minimizes the expression 

/o dz z 2A ~ 2 {(l - z 3 )[F'{z)\ 2 + A 2 z[F(z)} 2 } 



A 2 = 



To estimate it, use the trial function 



For A = 1 we obtain 



A 2 



lidz z^- 2 j^[F(z)] 2 
F = F a (z) = l-az 2 
6 - 6a + 10a 2 



a 2V37T - 6 In 3 + 4( % /3tt + 3 In 3 - 9) a + (12 In 3 - 13)a 2 
which attains its minimum at a ~ 0.24. We obtain 

A 2 sa Aq 24 ~ 1-27 

to be compared with the exact value A 2 = 1.245. The critical temperature is 



4tt 



+c 



4tt V A 



so for A = 1, T c « 0.225^/p, in very good agreement with the exact T c — 0.226^/p 
Similarly, for A = 2 we obtain 



1-fa+fa 2 



J 3-ln3 - ^ + (f -41n3)a+ (^ - ^ +ln3)a 2 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



whose minimum is A 2 « 17.3 (at a w 0.6) to be compared with the exact value A 2 = 16.754. The critical temperature 
in this case is T c w 0.117^/p, in very good agreement with the exact T c = 0.118^/p 5]. In fig. Q] we compare the 
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FIG. 1: The critical temperature in units of y/p vs. the dimension of the condensate A (eq. (|20l) ). The solid line is the analytic 
estimate (eqs. ()16[) and (JTTJ) ) and the almost indistinguishable dashed line is found by solving the field eqs. Q numericallly. 



analytic estimate of the critical temperature obtained from eqs. (fTI>)) and (jTTJ) with exact numerical results. The 
agreement between the two is excellent. 

Away from (but close to) the critical temperature, the field equation (|3|) for $ becomes 



$" = - V 4> (22) 

2rl A 1 - z 3 v 1 



where the parameter (Oa) /(2r3: ) is small. We may expand in the small parameter, 



\2 



r = A(i-«) + SxW + ... vm 

r + 2ri 



We deduce for the correction \ near the critical temperature 



with X (l) = x'(l) - 0. 

To find the temperature, we need 



From eq. (©, we deduce the ratio 



„2(A-l)p2/ \ 

1 + z + z z 



;1 2(A-l) F 2fl 

xi(0) = AC , C= / dz \> (25) 

1 + z + z A 



£=M 1 + ^ + -) (26) 



4 = a(i + ^ 

therefore the condensate near the critical temperature is 



T\ 2 /4tt x A 



Using the trial functions (fl7|) , for A = 1, we obtain from (|25[) . C ~ 0.54 and 7 « 11.4 to be compared with the exact 
7 = 9.3 @. Similarly, for A = 2, we find C « 0.07 and 7 w 133 to be compared with the exact 7 = 144 |5j. In fig. [5] 
we plot the analytic prediction (|27|) for the parameter 7 as a function of the dimension of the condensate A. Notice 
that 7 diverges as A -> 3. 
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FIG. 2: The parameter 7 that determines the condensate near the critical temperature vs. the dimension of the condensate A 
(eq. (2ZK). 

IV. LOW TEMPERATURES 

Turning to low temperatures, as T — > 0, we expect a simple scaling, ^ = ^>(bz), $ = where 6 — > 00. Then 

scaling z — » z/6 and letting 6 — > 00, the field equations ([3]) simplify, since the dominant contribution comes from the 
neighborhood of the boundary {z = 0) . Thus at low temperature we obtain the simplified system of equations 



z r\ 



,2A 



F 2 Z 2(A-1) $ = q 



(28) 



where we restored the original coordinate z (before scaling) . 

This system of coupled non-linear equations is to be solved subject to the boundary condition at the horizon 



3F'(1) + A 2 F(1) = 



(29) 



which is obtained from ([3]), © and (fT3")) . as well as those at the boundary, F(0) — 1 and F'(Q) — 0. 

For A = A_ < 3/2, numerical results indicate that F — > 1 as T — » 0. The solution to the field equation (|28p for $ 

is 



<S>(z)=Ar + Vz'K^(b A z A ) 



, A _ (OA) 

" Arf 



(30) 



The other solution is rejected because it is large at the horizon, contradicting the boundary condition ©. Notice 
that at the horizon $(1) ~ e~ b , which is an exponentially small error in the T — > (b — > 00) limit. 
From the behavior of $ (eq. (00])) near the boundary (z — 0), using ((6|, we deduce 



r(-ak) 



2 1+ 23 



(31) 



The field equation (|2"51) for F becomes 



F" + — — F' - A 2 bz(K^_ (b A z A )) 2 F = 



(32) 



Using perturbation theory, we obtain the solution 



F{z) = 1-A 2 b / dz 1 z' 
Jo 



dz" z"(K^(b A z" A )) 2 



(33) 
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Applying the boundary condition (|29|) at the horizon, we obtain 

A 2£2(2-A) 



For large 6, the denominator scales like fr 3_2A , showing that A ~ Then eq. (f5T|) implies p/r+ ~ 6 3 / 2 . Since 
(O a ) - (?- + 6) A (see eq. ©), it follows that 



Tc \T t , 

showing that the condensate diverges as (Oa) ~ T _A / 3 for A in the range 1/2 < A < 3/2 signaling the breakdown 
of the probe-limit approximation at low temperatures. 

Turning to the case A = A + > 3/2, notice that as we switch from A = A_ to A = A + , the boundary conditions 
at z = change, but not at the horizon. Thus, for a given to, the electrostatic potential $ has the same asymptotic 
behavior for both A + and A_. In terms of the scalar field, this implies F«l near the boundary (z = 0), whereas 



F(z) « - (36) 



a 

asymptotically (z > 1/6, where b 3> 1 is to be determined). Then in the asymptotic regime, eq. (|28p for $ has solution 
(c/. eq. (EOD) 

<S>(z)=Ar + WzK^(b 3 - A z 3 - A ) , b A = { ° A)a **~* (37) 
v y 2 ( 3~A) v ; (3 _ A)r A 

Notice that a singularity appears to develop at A = 3 indicating the onset of a quantum phase transition. 
Near the boundary, we deduce the estimate for the ratio p/r+, 

4 = - T -^^Ab (38) 

However, this is not a good estimate. In fact, it diverges at A = 5/2. We shall improve on this estimate by better 
accounting for the behavior of F near the boundary (where eq. (|36|) ought to be replaced by F w 1). 
The field equation (H for F is 

_ p „ 2(1-A) F , _ i (6 3-A z 3-A ))2f = (3g) 

z 2(3-A) 

Unlike with A < 3/2, in this case the term involving the electrostatic potential $ cannot be treated as a perturbation. 
By rescaling z — > z/b, eq. (|51?1) becomes 

-F" + 2{1 ~ A) F' -A 2 z(K i (z 3 ~ A )) 2 F = , A=4 (40) 

Z 2(3-A) 5 

With 

A~b (41) 

eq. (|40[) is independent of temperature and provides a good approximation to F and the corresponding eigenvalue A 
at T — 0. It ought to be solved in the interval (0, oo) subject to the boundary conditions F(Q) — 1, _F'(0) = 0, F 
as z — > oo (see eq. ([36]) ). These conditions determine the eigenvalue A. 
To estimate A, note that 

A 2 = ^ - — (421 

^dzz^[K_^_{z^)F{zW 

The cigenfunction F(z) minimizes this expression. We may substitute the trial function 

(ry\ 2A — 3 / 7 \ 2A — 3 

-) tanh(-) (43) 



1.8 2.0 2.2 2.4 2.6 2.8 

FIG. 3: The eigenvalue A denned in (|40[) as a function of A [solid line] compared with the estimate (142[) [dashed line]. 
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FIG. 4: The parameter a that determines the approximation (|43|l to the scalar field as a function of A. 



which obeys the correct boundary conditions. It interpolates smoothly between a constant value (F — 1) near the 
boundary and the power behavior (|36l) away from the boundary. The parameter a is fixed by minimizing the ratio 
(|4"2")l . Similar functions have been considered before [23, 24] but without a determination of the parameter a. In fig. [3] 
we compare the analytic estimate of the eigenvalue A using the trial functions (|43|) with exact numerical results. The 
agreement between the two is excellent. Notice that the limit A — > 3 is singular. A — > in this limit. Also, the 
parameter a labeling the trial function that minimizes (I42p . which is plotted in fig. 3] as a function of the dimension 
A, diverges as A — > 3. 

For A = 2 the minimum is obtained for a ps 0.8 which yields A ~ 1.92. In this case, A can be found exactly 
because eq. (|40)) can be solved analytically. We find explicitly 



F(z) 



2z 



which obeys the correct boundary conditions at z = 0. Demanding i 7, — > as z — ¥ oo then yields 

Jo 

showing that 



(44) 



(45) 



(46) 



where £o is the first root of the Bessel function Jq. Numerically, A = 1.91, in good agreement with our earlier estimate. 
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Restoring b, 



F(z) = —Y (Z )J (Z e 
so F(z) = 1 + 0(z 2 ), as desired. Away from the boundary (z > 1/6), 



-bz\ 



F(z) 



2bz 



which upon comparison with (|36l) yields 



7ry (&) 



0.8 



(47) 



(48) 



(49) 



in excellent agreement with our earlier estimate. 

To calculate the condensate at T = 0, we need a better estimate of the ratio (j3"5)l . To this end, we shall solve 
eq. (|28l) for $ perturbatively using (|37l) with A = 2 as our zeroth-order solution. We obtain 



We deduce the ratio 



a na 



„2fcz 



sin na 



: F{-a,l-l-a--e 2bz/a ) 



na z a 1 1 1 — a 

<* + -. + ^[ 

siri7ra 2 V 2 



2(1 - a) 



3 2b2/Q F(l-a,l,2-a;-f 



2bz/a 

(50) 
(51) 



=7ry (lo)/2 



improving the zeroth-order result ^fp/r + rj y^o b (eqs. (|3"5)) with A = 2 and (J3S])). Numerically, 5 — 0.58 and the 
condensate at T = is 



VW2 



1 + 



O.llSv 7 ^) 



7.9 



(52) 



in good agreement with the exact result (C'2) 1 ^ 2 = 8.3T C . 

To generalize to arbitrary A > 3/2, substitute the approximation to the function F a (z) (eq. (|4"3l ). 



F a (z) 



1 , z < a 

(a/z) 2A ~ 3 , z>a 



(53) 



into the field equation for the electrostatic potential $ and solve it to find a better approximation for $ than ([37 
After rescaling, z — > z/b, for z > a we obtain 



$ = $>(z) = .i&r+Vlfi: 1 (z 3 - A ) 

v ' ~ v 2(3-A) v ' 



whereas for z < a, 

$ = $<(z) = B+y/zI 1 
providing the estimate for p improving on (|38 



3- A 



2 



v 2A-3 



■B-y/zI_ i 



3- A 



,A 



v 2A-3 



3- A 



(54) 



(55) 



(56) 



P T(l + ^) V2Aa 2A " 3 / 

The coefficients B± are found by matching the two expressions at z = a. For A = 2 we obtain B+ ss — 1.736r+, 
B- s» 1.966r + , therefore p ps 1.43& 2 r 2 and ^/ (C 2 ) /7c 7.9, as before. Fig.[5]shows the dependence of the condensate 
on the dimension A. Our analytic estimate is in good agreement with the exact numerical value for most of the range 
of A. It becomes increasingly unreliable as A —> 3. This is expected from our estimate of the electrostatic potential 
(|55|) which is singular in the limit A — > 3. 
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FIG. 5: The condensate at zero temperature as a function of A. The solid line is our anaytic estimate whereas the dashed line 
is the exact numerical result. 



V. CONDUCTIVITY 



The conductivity on the boundary is found by applying a sinusoidal electromagnetic perturbation in the bulk of 
frequency ui obeying the wave equation 



d 2 A 
dr 2 



VA = u?A , V = 2/* 2 



(57) 



where A is any component of the perturbing electromagnetic potential along the boundary. Eq. (|5T[) is to be solved 
subject to ingoing boundary condition at the horizon 



A ~ e~^ r * - (1 - z)-W(3r+) 
as z — > 1 (r* — > — oo), where r* is the tortoise coordinate 



dr 



1 



f(r) 6r+ 



(1 - z) 3 

In- £- - 2V3tan 

1 — z A 



2 + z 



(58) 



(59) 



with the integration constant chosen so that the boundary is at =0. We shall solve this equation in the entire 
frequency spectrum. 

To this end, we shall replace the potential V with its average (V) in a self-consistent manner. We readily obtain 
the solution 



The other solution is discarded because it contradicts the boundary condition at the horizon 
conductivity 



The average value of the potential is found from 



J^VjA^ 
fLdr*\A(.r*)\ 2 



(60) 

We deduce the 
(61) 

(62) 



The integrals are well-defined if w has an imaginary part (which should be set equal to zero at the end of the 
calculation) . 
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For A < 3/2, in the low temperature limit the potential simplies to 

(Oa) 2 



V 



2(A- 



(63) 



where we used eq. (fT3"]) with F(z) s» 1. Moreover, since r + — > 0, the main contribution to the integrals in (R)2"j) is from 
the vicinity of the boundary where w —z/r + . We deduce the leading contribution 



(0a) 2 / °°dzz 2 ( A - 1 )|A(-z/r + )|- 
r 2(A-i) /-|A(-z/r + )P 



= T(2A- 1)(Oa) 2 -2t^oj 2 -(V) 



2(1- A) 



(64) 



which determines (V) implicitly as a function of to. We obtain the low-temperature high-frequency (oj > (Oa) 1 ^) 
conductivity 



T(2A- 1)(C A ) 2 (-2 Z )2( 1 - A ) 



2 A 



whereas for low frequencies, we have 



[2 2 (!- A )r(2A- 1)(C A ) 2 ] 1/A 



(65) 



(66) 



In particular for A = 1 cq. (j64|) can be solved for all frequencies and the expression (|65|) for the conductivity, which 
coincides with (1661) . is valid in the entire spectrum, 



(Oi 



a(u>) = \ 1 

V ijj- 

This expression is in excellent agreement with numerical results even down to low frequencies (uj <C (Oi))- 
For A > 3/2, the potential can be approximated by 



(67) 



Va(3- A) 2 6 2 r 2 (1 - z 3 )(bz) 2 ^-^ tanh 2 H^j 



2A-3 



(68) 



where b is given in ()37|1 . This expression can be used, as before, to find an estimate for (V). In particular, for A = 2, 
we obtain 



V = 1 + 2a\J V - lo 2 + 2a 2 (V - lu 2 ) 



if, 



where 



V 



(V) 
(0 2 )a 



(0 2 )c 



At high frequencies, this implies 



ct(w)rm/1 + 



(g 2 ) 2 
2w 4 



showing that a > 1 for to > y/ (O2), whereas as ui — > 0, V « 0.65, and the low-frequency conductivity is 

a(u) w 0. 



(69) 



(70) 



(71) 



(72) 



We shall improve on this estimate later by using a more accurate analytic technique which is better suited for low 
frequencies. We shall also obtain an exponentially small real part of the conductivity which survives in the limit 

At intermediate frequencies, we may expand around V = Co 2 = 1. We obtain 



<t{lj) 



a(0 2 



(73) 



11 



1.2 
1.0 
0.8 

Re(cr) 0.6 
0.4 
0.2 



0.0 





i 





0.0 



- Im(cr) 



0.5 1.0 1.5 2.0 2.5 




0.0 0.5 1.0 1.5 2.0 



V<0 2 > V<0 2 > 
FIG. 6: The real and imaginary parts of the conductivity at low temperatures for A = 2. 



for u/^f{Ch) w ^/a w 0.9. We may also use perturbation theory to go beyond the leading order. Treating 5V 
V — (02)01 as a perturbation, we obtain the wavefunction 



l + ^-^e-^M^- + ^(^l ; i + ^ ; 
,,2 



-2^/{0 2 )/a r* 



where 



/3 = iot\ 



(O 2 )o 



- 1 



We deduce the conductivity 



V ur cj 



(74) 



(75) 



(76) 



improving on the leading order expression (|73j) . The conductivity resulting from our analytic procedure (both real 
and imaginary parts) is plotted in fig. [5] for the entire spectrum. 

At low frequencies, the above expressions do not properly account for the boundary condition at the horizon. To 
this end, define 



A = (1 - z)- lu / {3r ^g(z) 

where Q is regular at the horizon (z = 1). The wave equation (|57[) reads 



3(1 - z A )Q" + (9z 2 - 2(1 + z + z 2 ) — )Q' - 



6* 2 (i + 2z) iUJ ( 2 + z )( 4 + z + z2 ) w2 
z 2 r+ 3(1 + z + z 2 ) r. 



(77) 



G = (78) 



At the horizon we may expand Q in a Taylor series. We deduce the boundary condition 

2* 2 (1) % 



3 — 2iuj/r + 3r + 



0(1) =0 



(79) 



At low temperature, for A = 1, we have <J/ ~ z. We may rescale z — ¥ z/b, where b — (C 1 )/r + and then let 
b — > 00. We obtain the approximate solution 



A = 



(80) 



12 



which is valid for low frequencies (w <C (Ci))- We deduce the conductivity 

'(^"•l+lEE — 

The ratio c+/c_ is found by applying the boundary condition (|79|) . We obtain 



c+ _ a - 3 _ 2a 2(2a 2 - 3) 



-2« 



c_ a + 3 

We deduce the low frequency expansion 

i(Oi) 



a(a + 3) 2 r + 



— + 0(u 2 ) , a 



(Oi) 

T+ 



cr(w) 



1 - 2 



o-3 2a _ 2(2a 2 -3) ^_ 2a 
'a + 3 



a(a + 3) 2 r + 



-+0{u 2 ) 



in agreement with the leading order result (joT)) . The DC conductivity is 

a _3<d) 



5Ra(0) ~ e- 2Q = e" A »/ T 



2tt 



For A = 2, at low temperatures eq. (|78| reads 



3^" g' 



3b 2 tanh 2 



bz 8uj 2 



0.48(d) 



G = 



whose general solution is given in terms of Legendre functions, 



(81) 



(82) 



(83) 



(84) 



(85) 



1 — tanh 



1 + tanh ^ 
We deduce the conductivity 



bz\ 



c+i 



P+ Q . - t - tanh — + c_P~ Q 



a J 



t(-l + Vl+4a 2 ) 



tanh — 
a J 



a(oj) s» i 



w 0.85 - 0.30^ 



(86) 



(87) 



The ratio c+/c_ is found from the boundary condition at the horizon (|79|) . At z w 1 we have tanh — w 1, so we may 
approximate 



i(-l + V , l+4Q 2 ) 



(tanh — ) = 



2±a/2 



a J r(l=pa) 



1 — tanh — 



We obtain 



5(1) 



.r(i-a) r(i + a) 

Applying the boundary condition (|79|) . we obtain 



c+ = c _2b r ( 1 ~ oQ 
c_ r(l + a) 



r(i - a) 



4iw 



c_ 6 



r(l + a) 



+ / e- b 



6r+ 



e 3r + (89) 



(90) 



showing that at low frequencies, 



-A s /T 



3vM©2 

2tt 



0.43^(d) 



(91) 



to be compared with our earlier estimate Q720 of the imaginary part which was obtained via a different, less accurate, 
analytic method. The real part is exponentially small and was not detected earlier. 

The above method can also be applied to other values of the dimension A if one replaces the potential by its 
self-consistent average (V). Then by solving the wave equation (|57l) using perturbation theory, we obtain the real 
part of the conductivity in the limit w — > and therefore the gap A g (3t<r(0) ~ e~ Afl / T ) for all values of the dimension 
A. 
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VI. CONCLUSION 



We have discussed analytic calculations involving holographic superconductors in the probe limit [5[ . These systems 
are labeled by the dimension A of the operator that condences below a certain critical temperature T c . We found 
approximate explicit solutions of the non-linear field equations in the bulk near the critical temperature as well as in 
the zero temperature limit. We obtained an analytic expression for the critical temperature in terms of an eigenvalue 
associated with the field equation of the scalar and showed that it was in good agreement with numerical results. At 
low temperatures, we showed that the condensate diverges as (0a) ~ T~ A / 3 for A < 3/2 signaling the breakdown of 
the probe approximation. For A > 3/2, we obtained an expression for the condensate at zero temperature in terms 
of an eigenvalue associated with the field equation of the scalar and demonstrated agreement with numerical results. 
Our method becomes unreliable in the limit A — > 3. We presented evidence that this limit is singular signaling the 
onset of a phase transition [l6| . We also calculated the conductivity analytically for various values of A and obtained 
good agreement with numerical results. In the DC limit we showed that the real part of the conductivity behaves as 
e -A s /T anc ] found analytic estimates of the gap A g . 

It would be interesting to extend our results beyond the probe limit by including back reaction to the bulk metric. 
Studying the resulting field equations will enable us to take the zero temperature limit without the obstruction of 
numerical instabilities. This will elucidate the nature of the ground state. Work in this direction is in progress. 
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